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Abstract 

We study the characteristics of the low mass pp enhancements near threshold in 
the three-body decays B + — > ppK + and B + — > ppir + . We observe that the proton 
polar angle distributions in the pp helicity frame in the two decays have the opposite 
polarity, and measure the forward-backward asymmetries as a function of the pp 
mass for the ppK + mode. We also search for the intermediate two-body decays, 
B + ~ > pA ++ and B + — > pA°, and set upper limits on their branching fractions. 
These results are obtained from a 414 fb -1 data sample that contains 449 xlO 6 BB 
events collected near the T(45) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider. 

PACS: 13.25.Hw 
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After the first observation of the charmless baryonic B meson decay, B + — > 
ppK + [1,2], many three-body baryonic decays were found [3,4,5,6]. The dom- 
inant contributions for these decays are presumably via the b — > s penguin 
diagram, shown in Fig. 1 (a); however B + — > ppn + is believed to proceed 
via the b — > u tree diagram as shown in Fig. 1 (b). One interesting feature 
of these decays is that the dibaryon mass spectra show enhancements near 
threshold as conjectured in Ref. [7]. Many theoretical explanations [8] have 
been proposed to describe these enhancements in the dibaryon system, which 
seem to be a universal feature of all charmless baryonic B decays. Study of 
the proton polar angular distribution for the dibaryon system in the ppK + 
mode [9] indicates a violation of the b — > s short distance picture [10]. Explicit 
predictions for the dibaryon mass spectra [11] and the angular distributions 
[12,13] for B + — > ppK + /n + became available after the experimental findings 
were reported. 

In this paper, we study the three-body charmless baryonic B meson decays 
B + — > ppK + and B + — > ppir + . The differential branching fractions as a func- 
tion of the dibaryon mass and the polar angle distributions of the proton in 
the dibaryon system are presented. We also search for intermediate two-body 
decays in ppir + three-body final states. This is motivated by the observations 
of two-body decays of charmed baryons [14]. Many predictions based on QCD 
sum rules [15], pole models [16] and a topological approach [17] indicate that 
B + — > pA ++ and B + — > pA° should be observable in the large data samples 
accumulated at the B-factories. 

We use a 414 fb" 1 data sample, corresponding to 449 xlO 6 BB pairs, collected 
with the Belle detector at the KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) 
collider [18]. The Belle detector is a large-solid-angle magnetic spectrome- 
ter that consists of a silicon vertex detector (SVD), a 50-layer central drift 
chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), 
a barrel-like arrangement of time-of-flight scintillation counters (TOF), and 
an electromagnetic calorimeter composed of CsI(Tl) crystals located inside a 
super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to detect K\ mesons 
and to identify muons. The detector is described in detail elsewhere [19] . Two 
inner detector configurations were used. A 2.0 cm beampipe and a 3-layer 
silicon vertex detector were used for the first sample of 152 xlO 6 BB pairs, 
while a 1.5 cm beampipe, a 4- layer silicon detector and a small-cell inner drift 
chamber were used to record the remaining 297 xl0 6 -B-B pairs [20]. 



The event selection criteria are based on information obtained from the track- 
ing system (SVD and CDC) and the hadron identification system (CDC, ACC, 
and TOF). All primary charged tracks are required to satisfy track quality cri- 
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Fig. 1. The possible leading (a) b — > s penguin diagram and (b) b — > u tree diagram 
for B + — > ppK + and S + — > decays, respectively. 

teria based on the track impact parameters relative to the interaction point 
(IP). The deviations from the IP position are required to be within ±0.3 cm 
in the transverse (x-y) plane, and within ±3 cm in the z direction, where the 
z axis is opposite the positron beam direction. For each track, the likelihood 
values L p , L K , and that it is a proton, kaon, or pion, respectively, are deter- 
mined from the information provided by the hadron identification system. The 
track is identified as a proton if L p / (L p + Lx) > 0.6 and L v j (L p + L n ) > 0.6, 
or as a kaon if L K /{L K + L n ) > 0.6, or as a pion if L n / (L K + L n ) > 0.6. For 
particles with momenta at 2 GeV/c, the proton selection efficiency is about 
84% (88% for p and 80% for p) and the fake rate is about 10% for kaons and 
3% for pions; the kaon selection efficiency is about 85% and the pion to kaon 
fake rate is about 2%; the pion selection efficiency is about 88% and the kaon 
to pion fake rate is about 11%. 

Candidate B mesons are reconstructed in the B + — > ppK + and B + — > ppn + 
modes. We use two kinematic variables in the center of mass (CM) frame to 
identify the reconstructed B meson candidates: the beam energy constrained 
mass Mbc = \f E^ eam — p B , and the energy difference AE = E B — E hcam , where 
-^beam is the beam energy, and Pb and E B are the momentum and energy, 
respectively, of the reconstructed B meson. The candidate region is defined 
as 5.20 GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.1 GeV < AE < 0.3 GeV. The 
lower bound in AE for candidate events is chosen to exclude possible cross- 
feed background from the decays with additional pions to the search modes, 
e.g. B + -> ppK* + . From a GEANT [21] based Monte Carlo (MC) simulation, 
the signal peaks in a signal box defined by 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 
and \AE\ < 0.05 GeV, and there is no peaking background except cross-feed 
events between the ppK + and pp7[ + modes. 

The background in the fit region arises dominantly from the continuum e + e~ — > 
qq (q = u, d, s, c) process. We suppress the jet-like continuum background 
events relative to the more spherical BB signal events using a Fisher discrim- 
inant [22] that combines seven event shape variables, as described in Ref. [23]. 
Probability density functions (PDFs) for the Fisher discriminant and the co- 
sine of the angle between the B flight direction and the beam direction in the 
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T(4S) rest frame are combined to form the signal (background) likelihood C s 
(Cb). The signal PDFs are determined using signal MC simulation; the back- 
ground PDFs are obtained from the sideband data with M^ c < 5.26 GeV/c 2 . 
We require the likelihood ratio 1Z = C S /(C S + C b ) to be greater than 0.75 
and 0.85 for the ppK + and ppir + modes, respectively. These selection criteria 
are determined by optimization of n s j^fn s + rib, where n s and rib denote the 
expected numbers of signal and background events in the signal box, respec- 
tively. We use the branching fractions from our previous measurements [9,4] 
in the calculation of n s and use the number of sideband events to estimate rib- 
If there are multiple B candidates in a single event, we select the one with the 
best x 2 value from the vertex fit. The fractions of multiple B events are about 
8% and 10% for the ppK + and ppn + modes, respectively. 

We perform an unbinned extended likelihood fit that maximizes the likelihood 
function, 

e -(N a +N b ) N 

L = — — J] [N s P s (M hCi , AE t ) + N b P b (M hCi , AE t )_ , 

i=i 

to estimate the signal yield in the candidate region; here P s (Pb) denotes the 
signal (background) PDF, iV is the number of events in the fit, and N s and iVj, 
are fit parameters representing the number of signal and background events, 
respectively. 



For the signal PDF, we use a Gaussian function to represent the signal Mt, c 
and a double Gaussian for AE with parameters determined by MC simula- 
tion. We then modify these parameters to account for the discrepancies be- 
tween data and MC using the ppK + signal events ( M pP < 2.85 GeV/c 2 ). 
With this correction, we can gain about 5% more signal yield. The continuum 
background PDF is taken as the product of shapes in Mb c and AE, which 
are assumed to be uncorrected. We use the parameterization first used by the 
ARGUS collaboration [24], /(M bc ) oc M hc Vl - x 2 exp[-f(l - x 2 )], to model 
the Mb c background, with x given by Mbc/.E'beam an d £ as a fit parameter. 
The AE background shape is modeled by a normalized second order polyno- 
mial whose coefficients are fit parameters. Because the ppn + mode can contain 
non-negligible cross-feed events from the ppK + mode, we include the ppK + 
MC cross-feed shape in the fit for the determination of the ppn + yield. The 
cross-feed from ppir + to ppK + is negligible. Figure 2 illustrates the fits of the 
B yields in a proton-antiproton mass region below 2.85 GeV/c 2 , which we 
refer to as the threshold-mass-enhanced region. The fitted B yields are 632 
^28 an d 184 tig, for the ppK + and ppn + modes, respectively. 
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Fig. 2. Distributions of AE (with M bc > 5.27 GeV/c 2 ) and M hc (with \ AE\ < 0.05 
GeV), respectively, for (a) ppK + and (b) ppir + modes with proton-antiproton pair 
mass less than 2.85 GeV/c 2 . The solid curves, solid peaks, and dashed curves rep- 
resent the combined fit result, fitted signal and fitted background, respectively. The 
dot-dashed curve indicates the ppK + cross-feed background in the fit to the ppir + 
mode. 



Since there are two different detector configurations and the detection effi- 
ciency is dependent on M pp , we separate the data sample into two sets and 
determine the B yields in bins of M pp , where the signal PDF is assumed to be 
the same for all M pp bins. We generate corresponding M C samples in order to 
estimate the efficiencies properly. The partial branching fractions are obtained 
by correcting the fitted B yields for the mass dependent efficiencies for each 
data set; they agree well with each other for the two data sets and these results 
are then combined to obtain the final results. 



The differential branching fractions as a function of the proton-antiproton 
mass for both ppK + and ppir + modes are shown in Fig. 3, and the measured 
branching fractions for different M pp bins are listed in Table 1 and Table 2. 
Note that we have defined the charm veto: the regions 2.850 GeV/c 2 < M pp < 
3.128 GeV/c 2 and 3.315 GeV/c 2 < M pp < 3.735 GeV/c 2 are excluded to 
remove background from B decay modes containing an rj c , J/tp, ip', Xco, or Xci 
meson. These results supersede our previous measurements [4,9] with better 
accuracy. The width of the ppn + mode is narrower than that of the ppK + 
mode and agrees better with the theoretical expectation [11]. The error bars 
show the statistical uncertainties only. The listed yield is the sum from the 
fits for two different periods; the listed efficiency is an effective one obtained 
by combining the two different detector configurations. 
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Fig. 3. Differential branching fractions for (a) ppK + and (b) ppir + modes as a func- 
tion of proton-antiproton pair mass. The solid curves are theoretical predictions [11] 
that are scaled to the observed charmless branching fractions. The two shaded mass 
bins, 2.85 < M pP < 3.128 GeV/c 2 and 3.315 < M pP < 3.735 GeV/c 2 , are not 
counted in the charmless signal yields since they contain contributions from the 
intermediate resonances 7] c ,J/ip and ip',XcO,Xci mesons, respectively. 



Systematic uncertainties are determined using high-statistics control data sam- 
ples. For proton identification, we use aA^ piT~ sample, while for K/tt iden- 
tification we use a D* + — > D tt + , D° — ► K~n + sample. The average efficiency 
difference for hadron identification between data and MC has been corrected 
to obtain the final branching fraction measurements. The corrections are about 
9% and 14% for the ppK + and ppn + modes, respectively. The uncertainties 
associated with the hadron identification corrections are estimated to be 4.2% 
for two protons and 1% for one kaon/pion. Tracking uncertainty is determined 
with fully and partially reconstructed D* samples. It is about 1% per charged 
track. The 1Z continuum suppression uncertainty is estimated from control 
samples with similar final states, B + — > J/ipK + with J /if) — > The un- 

certainties for TZ selection are 2.5% and 4% for the ppK + and ppn + modes, 
respectively. A systematic uncertainty of 2% in the fit yield is determined by 
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varying the parameters of the signal and background PDFs. The MC statisti- 
cal uncertainty is less than 2%. The error on the number of BB pairs is 1.3%, 
where we assume that the branching fractions of T(4S) to neutral and charged 
BB pairs are equal. The systematic uncertainties for each decay channel are 
summarized in Table 3. We first sum the correlated errors linearly and then 
combine them with the uncorrelated ones in quadrature. The total systematic 
uncertainties are 6.5% and 7.4% for the ppK + and ppir + modes, respectively. 

Table 1 

The B yields from AE — M^ c fits for the B + — > ppK + data sample, detection 
efficiencies and branching fractions (£>) in different M p p regions. 



Mpp (G 


eV/c 2 ) 


Yield 


eff( 
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B (10 


-6 


) 
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-2.0 


ao -°-11.0 
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< 2.85 

Table 2 

The B yields from AE 

efficiencies and branching fractions (B) in different M pp regions. 



M bc fits for the B^ 



ppir + data sample, detection 



Mpp (GeV/c 2 ) 


Yield 


eff(%) 


B (KT 6 ) 


1.876 - 2.0 


5i.3li:| 


27 


6 


0.4118$ ± 0.03 


2.0-2.2 


83.il! 2 : 2 


26 


6 


0.70lg;J§ ± 0.05 


2.2 - 2.4 


32.7+fl 


24 


5 


0.30lg;^ ± 0.02 


2.4 - 2.6 


14.81* 


22 


9 


0.14l8:8e±0.01 


2.6 - 2.85 


l-6iy 


23 


2 


0.02I884 ± 0.00 


2.85 - 3.128 


17.5l£ 2 


22 


5 


o.i7i8:8 6 7 ±o.oi 


3.128-3.315 


-0 5+ 4 - 3 
u -°-3.0 


22 


3 


-0.0118:83 ±0.00 


3.315-3.735 
3.735 - 4.0 


0.61^ 

c rr+5.6 
°-°-4.3 


21 
20 


3 
6 


0.0118:85 ± 0.00 
0.0618:85 ± 0.00 


4.0-4.6 


°- Z -6.7 


20 


8 


0.0918:8? ± 0.01 


4.6 - 5.15 


'• D -4.7 


15 


8 


-o.nl8:8? ± 0.01 


All 


l-78±g;f| ± 0.13 


with charm veto 


1.60lg;^ ± 0.12 


< 2.85 


1.57l^| ±0.12 
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Fig. 4. Differential branching fractions vs. cos 6 P in the proton-antiproton pair sys- 
tem for (a) B + — ► ppK + and (b) i? + — ► pp7r + . The solid curve is the theoretical 
prediction [12]. 
Table 3 

Systematic uncertainties(%) in the branching fraction for each decay channel. 



Source 


ppK + 


ppir + 


Tracking 


3.1 


3.2 


Proton Identification 


4.2 


4.2 


K/-7T Identification 


1.0 


1.0 


Likelihood Ratio Selection 


2.5 


4.0 


MC statistical error 


1.4 


1.8 


Fitting 


2.0 


2.0 


Number of BB pairs 


1.3 


1.3 


Total 


6.5 


7.4 



We study the baryon angular distribution in the proton-antiproton helicity 
frame at M pP < 2.85 GeV/c 2 . The angle 9 P is denned as the angle between the 
baryon direction and the oppositely charged meson direction in the proton- 
antiproton pair rest frame, i.e. this angle is determined by p and K~ /tt~ , or 
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by p and K + /tt + . We use the same likelihood method to estimate the B yield 
in each 9 P bin. Again, the signal PDF is fixed and the background shape is 
allowed to vary. The cos 9 P distributions, shown in Fig. 4, for the ppK + and 
ppn + modes have opposite trends. This distribution for the ppn + mode does 
not match the theoretical prediction [12], which is based on an extrapolation 
of the ppK + data using the perturbative QCD framework. However, it does 
agree with the naive short distance picture for afi^ti weak decay. Particles 
directly associated with b decay are more energetic and the particle containing 
the spectator quark is generally less energetic. After boosting to the proton- 
antiproton rest frame, the fast moving anti-protons and 7T + 's are back-to-back 
most of the time. However, the b — > s gluon process for the ppK + case seems 
to completely disagree with this short distance picture. The baryon with the 
spectator quark moves faster in the B rest frame. The same phenomenon 
has been observed in B° — >■ pAir~ [25] decays. Another theoretical prediction 
proposes a long distance effect, namely pp rescattering through a hypothetical 
baryonium bound state [13], in order to explain the violation of the short 
distance picture for the ppK + mode. Since this long distance effect should 
also occur for the ppn + case, it seems that further theoretical investigations 
are needed to simultaneously explain the behavior of both the ppK + and ppir + 
modes. 

Because we have enough B + — > ppK + signal events in the threshold enhance- 
ment region, we separate this region into five sub-regions. Fig. 5(a)-(e) shows 
the efficiency corrected B yield as a function of cos 9 P for these five sub-regions. 
We define the angular asymmetry as A Sp = ^ + ~^~ , where N + and AL stand 
for the efficiency corrected B yield with cos 9 P > and cos 9 P < 0, respectively. 
The measured angular asymmetry as a function of M pp is shown in Fig. 5(f). 
It is interesting to see that there is a clear trend, which indicates that the 
relative contributions from two (or more) competing decay amplitudes are 
changing in this mass range. The measured average A 6p value of the threshold 
enhancement is given in Table 4. The systematic error, ~ 0.03, is determined 
by checking the B + — > J/ipK + (J/ip — > At + /i~) sample and the continuum 
background in B + — > ppK + where no asymmetry is expected. The observed 
A 9p 's are 0.02 ± 0.01 for B + -> J/^K + and 0.00 ± 0.02 for the continuum 
background. 



We also search for the intermediate two-body decays, B + — > pA ++ (A ++ — > 
pir + ) and B + — > pA° (A — > pn + ), from the ppn + three-body final state. 
Events with M pw < 1.4 GeV/c 2 are selected. No significant signals are found 
from the likelihood fit in those decay chains. We observe 59 and 86 events in 
the signal box; the expected numbers of background events from the fits are 
73.0 ± 1.6 and 81.4 ± 1.6 for B + -> pA ++ and B + -> pA°, respectively. We 
set upper limits on the branching fractions at the 90% confidence level using 



11 




o300 - 



LD200 - 



100 



(b) 



-1 -0.5 0.5 1 

Cose n 



&300 - 

CD 
"O 
<D 

£200 

0) 
> 

100 



(e) 



-1 -0.5 0.5 1 
Cose„ 



(c) 



100 



-1 -0.5 0.5 1 
Cose n 




2 2.25 2.5 2.75 

M - (GeV/c 2 ) 



Fig. 5. Efficiency corrected B yield vs. cos9 p for (a) M p p < 2.0 GeV/c 2 , (b) 
2.0 < M pP < 2.2 GeV/c 2 , (c) 2.2 < M P p < 2.4 GeV/c 2 , (d) 2.4 < M pP < 2.6 
GeV/c 2 , and (e) 2.6 < M P p < 2.85 GeV/c 2 ; (f) the measured angular asymmetries 
(A0 ) for these five mass regions near threshold. 

the methods described in Refs. [26,27] where the 7.4% systematic uncertainty 
for B + — > ppir + is taken into account. The results are B(B + — > pA ++ ) < 
0.14 x 1CT 6 and B(B + — > pA°) < 1.38 x 1CT 6 . These numbers are smaller than 
the theoretical expectations but agree with other experimental findings [28]. 

Since there is a prediction [29] that direct CP violation in B + — > J/ipK + is at 
the 1% level, it is quite possible that this effect could be magnified due to the 
interference [30] between the resonance and the threshold enhancement. We 
define the charge asymmetry A C p as (N b — A^)/(A^ + iVg) for the ppK + and 
ppn + modes, where Nb(Ni) stands for the efficiency corrected B~ (B + ) yield. 
The selection criteria for J ftp (rj c ) and related consistency checks have been 
reported in Ref. [31]. We adopt the same criteria and assume the signal PDFs 
are the same for both B~ and B + samples. The results from the likelihood 
fits are listed in Table 4 for various mass/resonance regions. No significant 
charge asymmetries are found. The systematic uncertainty is assigned using 
the measured charge asymmetry for sideband data and is found to be —0.01 ± 
0.01. 



In summary, using 449 x 10 6 BB events, we measure the mass and the angu- 
lar distributions of the proton-antiproton pair system near threshold for the 
ppK + and ppii + baryonic B decay modes. These results supersede our previous 
measurements [4,9] with better accuracy. The width of the threshold enhance- 
ment in the ppir + mode is narrower than that of the ppK + mode and agrees 
better with the theoretical expectation [11]. The proton polar angular distri- 
butions of the ppK + and ppn + modes have opposite trends. This shows that 
the b — > s and b — > u processes are kinematically different at short distance. 
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We also search for intermediate two-body decays in the ppn + final states; no 
significant signals are found. 

Table 4 

Summary of the results in the mass region M p p < 2.85 GeV/c 2 and in other in- 
termediate resonance regions. Y is the fitted signal yield (or upper limit at 90% 
confidence level), B is the branching fraction, Aq is the angular asymmetry and 
Acp is the charge asymmetry. 



Mode 


Y 


B (KT 6 ) 


A e 


A C p 


ppK + 


63211 


5.001^2 ± 0.32 


0.45 ± 0.05 ± 0.03 


-0.02 ±0.05 ±0.02 


VcK+ 


158111 






-0.16 ±0.08 ±0.02 


J/iJ;K + 


236+ljj 






0.09 ±0.07 ±0.02 


PpTT + 


184i» 


1.57l^±0.12 


-0.47 ±0.12 ±0.03 


-0.17 ±0.10 ±0.02 


pA ++ 


< 7.5 


< 0.14 








< 25.9 


< 1.38 
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